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ABSTRACT 

We examine the hypothesis that hydrostatic pressure alone determines the ratio of atomic to molec- 
ular gas averaged over a particular radius in disk galaxies. The hypothesis implies that the transition 
radius, the location where the ratio is unity, should always occur at the same value of stellar surface 
density in all galaxies. We examine data for 28 galaxies and find that the stellar surface density at 
the transition radius is indeed constant to 40% at a value of 120 Mq pc~^. If the hypothesis can 
be confirmed at all radii within a large range of galaxy types and metallicities, combining it with 
the observed constancy of the star formation rate with H2 surface density may enable a physically 
motivated star formation prescription with wide applicability. 
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1. INTRODUCTION 

Current ideas about the formation of GMCs present 
a nagging puzzle. Some authors have suggested that 
GMCs form by agglomeration of pre-existing molecular 
clouds (e.g. Scoville & Hersh 1979; Kwan & Valdes 1983). 
Others have argued that GMCs form primarily from 
atomic gas through some sort of instability or large-scale 
shock (e.g. Woodward 1976; Blitz & Shu 1980; Engar- 
giola, et al. 2003). In principle, both modes may occur: 
in galaxies that are predominantly atomic, GMCs might 
form by one process, and in galaxies that are primarily 
molecular, the GMCs might form by another. If so, two 
different CMC formation mechanisms could be at work 
in the same galaxy because many galaxies show a transi- 
tion from being predominantly molecular at their centers 
to being predominantly atomic in their outer parts (e.g. 
Wong & Bhtz 2002; Heifer et al. 2003). For example, in 
the outer parts of spiral galaxies there is so little molec- 
ular gas (e.g. Dame et al. 1987; Dame 1993; Heyer et 
al. 1998) that making GMCs from atomic gas seems to 
be the only available formation pathway. On the other 
hand, the centers of many galaxies are so overwhelm- 
ingly molecular (e.g. Mauersberger et al. 1989; Young 
et al. 1995; Wong & Blitz 2002, Heifer et al. 2003) that 
it is implausible for the inner galaxy molecular clouds 
to form from anything other than pre-existing molecular 
gas. The stars in these galaxies form only from molecular 
gas, yet there is no obvious change in the star forming 
properties across the molecular /atomic transition region. 
How could it be that molecular clouds that form by two 
independent processes would show no obvious difference 
in their star formation properties? 

One possibility is that the process of cloud formation is 
independent of whether the preexisting gas is atomic or 
molecular. That is, the ratio of atomic/molecular gas de- 
pends on some other factor, such as the ambient pressure. 
Because the rate of star formation d epends only on the 
amount of molecular gas l)Wong & Blit z 2002') . whatever 
determines the molecular gas surface density, determines 
the variation of star formation within a galaxy. 
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In this paper, we consider the hypothesis that the 
molecular gas fraction in a galaxy disk is determined by 
the mean hydrostatic pressure at a particular radius. We 
show that if the hydrostatic pressure is the only parame- 
ter determining the molecular gas fraction, one predicts 
that the radius at which the atomic and molecular gas 
surface densities are the same, the transition radius, oc- 
curs at a constant value of the stellar surface density, . 
Remarkably, in a sample of 28 galaxies, is found to 
be constant to about 40% at the transition radius, even 
though the observed variation of S* in these galaxies is 
at least 3 orders of magnitude. 

2. BACKGROUND 

Several authors have previously suggested that gas 
pressure determines the molecular fraction at a given ra- 
dius in a galaxy. Spergel & Blitz (1992), for example, ar- 
gued that the extraordinarily large molecular gas fraction 
at the center of the Milky Way is plausibly the result of 
the very high hydrostatic pressure in the Galactic bulge. 
Elmegreen (1993) suggested on theoretical grounds that 
the ratio of atomic to molecular gas in galactic disks 
results from both the ambient hydrostatic pressure as 
well as the mean radiation field. The dependence on 
the pressure is steeper than that of the radiation density 
{fmoi oc p2.2j-i.i^ g^j^i^ ought to bc the dominant factor. 
Observationally, Wong and Blitz (2002) showed that the 
radial dependence of the atomic to molecular gas ratio in 
seven molecule rich galactic disks can be understood to 
be the result of the variation in interstellar hydrostatic 
pressure (with fmoi cx P°-*). 

Let us assume, therefore, that A^(H2)/A^(H I), the ra- 
tio of H2 column density to H I column density, is deter- 
mined only by the midplane hydrostatic pressure, Pext- 
In an infinite disk with isothermal stellar and gas lay- 
ers, and where the gas scale height is much less than the 
stellar scale height, as is typical in disk galaxies, to first 
order: 
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where Sg is the total surface density of the gas, Vg is the 
velocity dispersion of the gas, is the midplane surface 
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density of stars, and pg is the midplane surface density of 
gas. The first term on the right is due to the hydrostatic 
pressure of the gas in the steUar potential; the second 
term is due to the self-gravity of the gas. 

In most galaxy disks, is much larger than pg when 
averaged over azimuth, except in the far outer parts 
of a galaxy where the stars become quite rare. In the 
solar vicinity, for example, — 0.1 Mq pc~^ (e.g. 
Binney & Merrifield 1998), but pg ~ 0.02 M© pc'^ 
(e.g. Dame 1993). For a self-gravitating stellar disk, 

S* = 2\/2/5*ft.*, where ft,* is the stellar scale height and 
ft* = (w*^/47rGp,)*'-^. Thus, neglecting pg, Equation 1 
becomes: 

Pe.* = 0.84(GI]*)°-^S3^£2_ (2) 

where Sg is the total surface density of the gas, Vg is the 
velocity dispersion of the gas, p* is the midplane surface 
density of stars, and pg is the midplane surface density of 
gas. Direct solution of the fluid equations by numerical 
integration shows that this approximation is accurate to 
10% for S* > 20 Mq pc^^ (where p* pg), which covers 
the range of stellar surface densities in this study. 

We choose to express the midplane pressure in the form 
of Equation [3 since there is good evidence that both ft* 
and Vg are constant with radius in disk galaxies. Further- 
more, because of the weak dependence of Pext on ft* , and 
the small variation of ft* measured among galactic disks, 
we expect variations of ft* to have little effect on Pext- 
The constancy of the stellar scale height within galaxies 
was demonstrated by Ivan der Kr^iit fc Searld l)1981albf) 
and has been c onfirmed in other edge-on galaxies (e.g. 
IFrv et al.llT999|) . While there is some evidence that the 
stel lar scale height flares at large radius in some galax- 
ies (jNaravan fc Jog 2002"!. this is only found only at the 
edges of stellar disks in regions where the stellar, gaseous 
and dark matter components of the disk make compara- 
ble contributions to the potential. Adopting a constant 
stellar sca le height i s furt her supported by the obser- 
vations of 'Bo ttemal ljl993D . who shows that the stellar 
velocity dispersion follows an exponential distribution 
with scale length twice that of the stellar surface den- 
sity in disks. This observation suggests that the stel- 
lar disk is distributed vertically in a sech^z profile with 
cr* = (\/27rGI]*ft*)i/2. 

While the stellar component of galactic disks can be 
approximated with a constant scale height, the gas com- 
ponent is better described by a constant velocity dis- 
persion, a„, whi ch is ob serve d in face-on ga laxies (e.g. 
iShostak fc van der Krifi^[l9 84: Dick ev et al2f l990i and 
the Milkv Wav l|Burto'nl ¥971: Malhotral Il99fi1 . The 
Milky Way observations show that the gas scale height 
decreases as the stellar surface mass density increases 
in a manner co nsistent the gas remaining isothermal 
l)Malhotral 1 99,'t() . This ensemble of observations suggests 
that ag = 8 km s~^ characterizes the H I velocity disper- 
sion in the stellar dominated regions of galactic disks. 

By assumption, 

7V(H2)/iV(H I) = fiPext) (3) 

Thus, since (ug/\/ft7) is approximately constant within 
galaxies, then 

7V(H2)/7V(HI)-/[Pe.t(Sg,I]*)] (4) 



The mass surface density of atomic gas, Shi, is rea- 
sonably constant across the inne r portions of galac- 
tic disk (e.g. iWevers et al.l Il986t iCavatte et aD 119941 
iWong fc 3111^120021) . While some changes of Sri with 
radius are observed, the variation is small compared 
to the changes observed in stellar surface mass density 
and molecular surface mass density. We may therefore 
adopt a single value for the surface density of atomic 
gas in a galaxy across the stellar disk. iV(H I) satu- 
rates at a value of about 1 x 10 cm" -^ ^Wong fc Blitd 
|2005|). Thus, N(H2)/N(H I) has a value of unity when 
2iV(H2)-l-iV(H I) has a value of - 2 x lO^^ cm'^. Equa- 
tion[21 therefore implies that the radius where the atomic 
and molecular surface densities are equal in spiral galax- 
ies depends only on E* if Ug/ft* is constant. 

3. RESULTS 

We use BIMA SONG fH elfer et alJl200^ . an interfer- 
ometric imaging survey of the CO in 44 nearby spirals, 
to determine the distribution of the molecular gas. To 
determine the H2 surface density, we use a conversion fac- 
tor of iV(H2)/TA(C0) Aw = 2x1020 cm-2(K km s^^)-^. 
For the H I, we adopted a single value that characterizes 
the column density across the galactic disk. These values 
were drawn from maps of galaxies reported in the liter- 
ature. Adopted values and the corresponding references 
appear in Table H] alo ng with the o r ientat ion parameters 
and distances used in iHelfer et alJ l)2003H . If no H I ob- 
servations have been reported for the galaxies, we used 
a value of 8 Mq pc"^, which corresponds to a surface 
density of 1 x 10^^ cm"^, typical for the stellar disks of 
galaxies. 



TABLE 1 

Adopted Galactic Parameters 



Galaxy 




Inc. 


P.A. 


^HI 




Reference 


Name 


(Mpc) 


(°) 


(°) 




(M^/pc2)3/2 




NGC 628 


7.3 


24 


25 


4.9 


52 


Kamohuis & Briees (1992) 


NGC 1068 


14.4 


33 


13 


13.5 


161 


Brinks et al. fl997l 


IC 342 


3.9 


31 


37 


4.0 


59 


Crosthwaite et al. ("20011 


NGC 2903 


6.3 


61 


17 


5.0 


63 


Wevers et al. fl986") 


NGC 3184 


8.7 


21 


135 


8.0 


81 




NGC 3351 


10.1 


40 


13 


5.0 


56 


Schneider (1989') 


NGC 3368 


11.2 


46 


5 


4.0 


46 


Warmels ^19881 


NGC 3521 


7.2 


58 


164 


8.0 


75 




NGC 3627 


11.1 


63 


176 


2.2 


23 


Zhane et al. ("1993") 


NGC 3726 


17.0 


46 


10 


7.0 


62 


Wevers et al. ('1986') 


NGC 3938 


17.0 


24 





7.5 


67 


van der Kruit & Shostak (19821 


NGC 4051 


17.0 


41 


135 


4.0 


36 


Wone & Blitz (20021 


NGC 4258 


8.1 


65 


176 


4.5 


53 


Wevers et al. (19861 


NGC 4303 


15.2 


27 





8.6 


96 


Cavatte et al. (19941 


NGC 4321 


16.1 


32 


154 


5.0 


49 


Wone & Blitz (20021 


NGC 4535 


16.0 


45 


28 


4.3 


38 


Cavatte et al. (19941 


NGC 4569 


16.8 


62 


23 


5.4 


51 


Cavatte et al. (19941 


NGC 4579 


16.8 


37 


95 


3.2 


27 


Cavatte et al. (19941 


NGC 4736 


4.3 


35 


100 


8.0 


123 


Wone & Blitz (20021 


NGC 4826 


4.1 


54 


111 


10.0 


118 


Braun et al. (19941 


NGC 5005 


21.3 


61 


65 


8.0 


116 




NGC 5033 


18.6 


62 


170 


8.0 


62 


Wone & Blitz (20021 


NGC 5055 


7.2 




105 


6.3 


72 


Wone & Blitz (20021 


NGC 5248 


22.7 


43 


110 


8.0 


93 




NGC 5247 


22.2 


29 


20 


8.0 


73 




NGC 5457 


7.4 


27 


40 


6.3 


81 


Wone &c Blitz (20021 


NGC 6946 


5.5 


54 


65 


7.8 


95 


Tacconi Youn? (19861 


NGC 7331 


15.1 


62 


172 


8.0 


87 





The stellar surface densities are determined using re- 
duced K-h&nd images from the 2MASS Large Galaxy 
Atlas l|Jarrett et alJ 12003(1 . We adopted a constant 
JC-band mass-to-light ratio of Ml/Lk — 0-5AIq/Lq 
l|Bell fc de Jonell2001(l . Our final sample consists of only 
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those galaxies with both 2MASS data and CO detec- 
tions. Of these, 22 galaxies have H I surface densities in 
the literature and 6 galaxies do not. The stellar and gas 
surface densities are corrected to face-on values using the 
orientation parameters listed in Table ^ 

We define the transition stellar surface density (S*,t) 
where iV(H2)=A^(H I), as the median value of the stellar 
surface density at all positions for which 0.9A^(H I) < 
^(H2) < l.liV(Hl). We associate this surface mass 
density with the transition radius Rt in the galaxy by 
finding the radius where the azimuthally averaged stel- 
lar surface mass density S*(i?gai) equals the transition 
surface density, In Figure ^ we plot the transition 

radius against the transition stellar surface density, using 
the adopted distances in Tabled We also list the values 
of Shi-v/STJT are listed in Table as a check since Pext 
depends directly on this quantity. 



10^ 



1 0-1 r 



o 



101 - 




« CO, HI 

o CO. No HI 
^Milky Way 
■Jk-M33 



: ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ : 

2 4 6 8 10 

Rtrans [kpc] 

Fig. 1. — Plot of the median stellar surface mass density where 
N(Ji2) ~ N{R I) as a function of where this surface density occurs 
in the galaxy. For galaxies with measured H I densities (filled 
circles, 22 galaxies), the range of stellar surface densities is plotted 
as an error bar, running between the values of the 1st percentile 
and the 99th percentile of surface density in the galaxy. Galaxies 
without H I measurements are plotted as open circles (6 galaxies). 
The transitions stellar surface density is remarkably constant and 
has a mean value of 120 it 10 Mq pc"'^ for the 22 galaxies with 
both CO and H I data. Points for the Milky Way (MW) and M33 
are also plotted. 



Figure^shows a remarkable constancy of S^^t expected 
if iV(H2)/A^(H I) is determined by pressure only; the 
mean value is 120 ± 10 M© pc~^. Also plotted in the 
figure is the range of measured for each individual 
galaxy. Formally, the dispersion in the mean value of 
E*^t, is only 40%, yet the the range of to which the 
2MASS survey is sensitive in these galaxies varies by al- 
most 3 orders of magnitude. The scatter in Eniy^E^ 
given in Table is only 60%. The range of galactocen- 
tric distance where the transition radius occurs varies by 
more than an order of magnitude. Apparently, the con- 
stancy of E» J is not due to a small range in the observed 
properties of the galaxies. The small scatter in t also 
suggests that various assumptions such as that of a con- 
stant value of ft.* both within and among galaxies are 



justified. 

As a check, one can calculate the value of E*^t for the 
Milky Way, scaling the measured E* at the distance of 
the Sun (35 Mg pc'^; Binney & Merrifield 1998), and a 
radial scale length for the stars of 3 kpc (Spergel, Malho- 
tra & Bhtz 1996; Dehnen & Binney 1998). The transition 
radius for the Milky Way occurs at a galactocentric dis- 
tance of about 4 kpc (Dame et al. 1993), or about 1.3 
scale lengths inward of the Sun. This converts to a stellar 
surface density of 132 M© pc~^, in good agreement with 
the determinations in other galaxies. The trend can also 
be checked in M33 using the data for fm.nl presented in 
iHever. Corbelli. Schneider, fc Youiiel l)2004D . which gives 
E,,t = 19OM0 pc-2 where Ehi = Eh2 {Rgai = 300 pc). 

There appears to be a small but significant decrease 
in the transition surface density with radius, which may 
be due to a breakdown of the assumptions of a constant 
mass-to-light ratio, Vg and /i* with radius, or that pres- 
sure alone determines iV(H2)/iV(H I). A linear fit to the 
data in Figure 1 gives logE* = (2.36 ± 0.05) - (0.06 ± 
0.01)(i?/l kpc)MQ pc~^ with errors given by the scatter 
in the data around the trend. 

4. DISCUSSION 

Figure 1 is consistent with the hypothesis that the 
mean hydrostatic pressure determines the ratio of atomic 
to molecular gas at a given radius in a disk galaxy. The 
small scatter in the mean value of E* suggests that glob- 
ally, the pressure may be the only important factor in de- 
termining the ratio of atomic to molecular gas. But the 
variation in the hydrostatic equilibrium of a disk is ex- 
pected to be rather smooth. Why then do galaxies show 
so much azimuthal variation in the molecular gas, and by 
implication in the atomic-to-molecular gas ratio (Heifer 
et al. 2003)? Significant variations in the interstellar 
pressure can result from a variety of causes such as spi- 
ral shocks and explosive events (e.g. supernovae). Thus 
large pressure variations can occur on all scales locally 
even if the mean hydrostatic pressure varies only slowly. 
Furthermore, even if the hydrostatic pressure drives the 
ratio globally, locally the radiation field can be impor- 
tant in determining ho w much of a m olecular cloud can 
remain neutral (HoUcn bach fc Tielens 1999). Thus, sig- 
nificant deviations from the mean molecular abundance 
can be produced by both pressure and radiation varia- 
tions. 

Pressure or density? The rate of formation of molecu- 
lar gas is thought to be dependent on the local gas density 
in the chemical reactions that produce H2 and CO. Are 
we then using pressure as a surrogate for the mean gas 
density in determining A^(H2)/A^(H I) as a function of 
galactocentric distance? Pressure, as defined in Equa- 
tions 1 and 2, is taken to be Pgfg where Gg includes 
both thermal and turbulent contributions. Because ag is 
measured to be typically 7-8 km s~^and the gas tem- 
perature of the cold gas layer where most of the gas mass 
resides is typically < 100 K, the thermal pressure is only a 
small fraction of the turbulent pressure. Observationally, 
(jg is constant with radius (see §2), thus Pext oc pg] vari- 
ations in pressure are effectively the same as variations 
in density. We choose to describe the functional depen- 
dence in terms of pressure rather than density because 
it is directly measurable through Equation 2 (assuming 
/i» is known) , whereas the density is inferred and not di- 
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rectly measurable on galactic scales. But it should be 
kept in mind that given the measured constancy of ag, 
we cannot distinguish the effects of pressure from those 
of density. 

What pressure is implied by S*.t? In the Milky Way, 
the value of /i, is about 300 pc (Binney & Merrifield 
1998), Vg is about 7 km s~^ (Dickey & Lockman 1990), 
Eg is 8.6 Mq pc~^ and is 132 Mq pc"^ Using Equa- 
tion |21 Pext/k = 1.5 xlO^ cm~^ K after correcting for 
helium. This value is still an order of magnitude below 
the mean internal Pint/k ~ 3 xlO^ cm~^ K for GMCs 
that have t ypical surface densities of ~ 100 Mq pc~^ 
(jmEiliin^). Because P„t > Pe^t, GMCs that survive 
for more than a crossing time, ^ 10^ y, must be self- 
gravitating. Using Equation 2 to scale Pext to the inner 
regions of disks suggests that GMCs are self-gravitating 
over nearly the entire disk. 

What are the implications for understanding star for- 
mation on galactic scales? If the global atomic-molecular 
transition is governed by pressure across a wide range 
of galaxies, it will be possible to develop a perscrip- 
tion for star formation on global scales that is physi- 
cally well motivated. Non-thermal radio-continuum is 
tight ly correlated with the far IR emission in galaxies 
(e.g. iCond on 1992), implying that the radio-continuum 
is a good extinction-fr ee indicator of the st ar formation 
rate in spiral galaxies. iMurgia et alJ ()2002D have shown 
that for a sample of 180 spiral galaxies, the ratio of radio- 
continuum to CO emission is constant to within a fac- 
tor of 3, suggesting that the star formation efficiency of 



molecular clouds averaged over galactic scales is constant 
at about 3.5%. Therefore, if the relationship between 
pressure and A^(H2)/iV(H I) (i.e the function / in Eqs. 
1 and 2) can be determined for all galaxies, or if the 
variation in / can be found for different galaxy types, 
then it will be possible to determine the star formation 
rate in galaxies by measuring the stellar and gas surface 
densities only. Furthermore, it will be possible to obtain 
reliably the star formation rate from simulations, since 
the turbulent gas pressure can be directly calculated. In 
addition, if the variation in / can be measured for galax- 
ies of low metallicity, then determining the star formation 
rate can be extended to high z. The measurement of / 
will be the subject of a future paper. 

We suggest, then, that GMCs can form from either 
preexisting atomic or molecular gas depending on the 
dominant state of the diffuse interstellar medium at a 
particular radius in a galactic disk. That dominant state 
is determined by the hydrostatic pressure, modified by 
local perturbations such as density waves, supernova re- 
manants, etc. 
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